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Guiding  Supersonic  Projectiles  using  Optically  Generated  Air  Density  Channels 

Phillip  Sprangle  and  Luke  A.  Johnson1 

Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington  DC  20375 
'University  of  Maryland,  College  Park,  Maryland 

Abstract 

We  investigate  the  feasibility  of  using  optically  generated  channels  of 
reduced  air  density  to  provide  trajectory  correction  (guiding)  for  a  supersonic 
projectile.  It  is  shown  that  the  projectile  experiences  a  force  perpendicular  to  its 
direction  of  motion  as  one  side  of  the  projectile  passes  through  a  channel  of  reduced 
air  density.  A  single  channel  of  reduced  air  density  can  be  generated  by  the  energy 
deposited  from  filamentation  of  an  intense  laser  pulse.  We  propose  changing  the 
laser  pulse  energy  from  shot-to-shot  to  build  longer  effective  channels.  Current 
femtosecond  lasers  systems  with  multi-millijoules  pulses  could  provide  trajectory 
correction  of  several  meters  on  5  km  trajectories  for  sub-kilogram  projectiles 
traveling  at  Mach  2  to  Mach  3. 

I.  Introduction 

The  purpose  of  this  paper  is  to  consider  the  possibility  of  guiding  supersonic  projectiles  by 
forming  channels  of  reduced  air  density  in  front  of  the  projectile  as  shown  in  Fig.  1.  A  reduced 
air  density  region  on  one  side  of  the  projectile  would  create  a  surface  pressure  imbalance,  resulting 
in  a  net  transverse  force  on  the  projectile  and  thus  the  possibility  of  trajectory  correction,  i.e., 
guiding.  Air  density  channels  are  formed  by  creating  plasma  filaments  ahead  of  the  projectile  with 
high-intensity  lasers  pulses.  The  plasma  filaments  consist  of  hot  electrons  which  collisionally  heat 
the  surrounding  air  forming  a  density  depression.  The  filaments  are  generated  by  a  dynamic 
balancing  between  the  laser’s  nonlinear  self-focusing  and  plasma  defocusing.  From  shot-to-shot, 
the  location  of  the  laser  focus  can  be  controlled  by  a  combination  of  adjusting  the  initial  conditions 
of  the  laser  pulse,  such  as  the  pulse  energy.  A  flying  focus,  consisting  of  successive  laser  pulses 
focusing  at  different  locations,  can  produce  air  density  channels  ahead  of  the  projectile  over  a 
distance  much  larger  than  a  single  filament. 

There  is  extensive  literature  on  using  plasmas  for  supersonic  flow  control,  specifically  for 
the  purpose  of  drag  reduction  (Touchard  2008).  In  these  systems,  plasmas  are  generated  as  either 
a  means  to  deposit  energy  and  disrupt  bow  shocks  (Kim  201 1)  or  as  a  means  to  modify  the  surface 
flow  and  avoid  laminar  flow  separation  (Touchard  2008).  Our  work  is  distinct  from  this 
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body  of  literature  because  it  is  focused  on  the  feasibility  of  projectile  guiding  instead  of  drag 
reduction.  The  concept  of  using  laser  energy  deposition  for  shock  modification  and  projectile 
guiding  was  previously  proposed  in  a  U.S.  patent  (Kremeyer  2003)  but  the  feasibility  and  physics 
of  the  mechanism  were  not  covered. 

The  paper  consists  of  three  sections  dealing  with  a)  the  degree  of  projectile  deflection  given 
by  an  air  density  channel,  b)  the  formation  and  characteristics  of  air  density  channels  and  finally 
c)  the  procedure  for  forming  a  flying  focus. 


FIG.  1.  Configuration  of  the  laser  pulse  system,  plasma  filament,  and  air  density  channel  with 
respect  to  the  supersonic  projectile. 

II.  Supersonic  Projectile  Deflection 
A.  Forces  on  Projectile 

To  determine  the  effect  of  an  air  density  channel  on  the  projectile’s  trajectory,  the  total 
force  on  the  projectile  must  be  calculated.  As  the  projectile  moves  through  the  air,  it  is  only 
interacting  with  the  gas  through  the  air  pressure  over  its  surfaces.  Therefore,  in  principle,  the  net 
force  is  found  by  integrating  the  pressure  over  the  projectile’s  surface.  When  the  net  force  is 
decomposed  into  components  anti-parallel  and  perpendicular  to  the  projectile’s  direction  of 
motion,  it  is  referred  to  as  the  drag  and  lift,  respectively.  The  force  perpendicular  to  the  projectile’s 
motion  will  be  referred  to  as  the  “transverse  force”  to  emphasize  that  its  purpose  is  for  deflecting 
the  projectile  and  not  solely  for  lift.  The  projectile’s  drag  force  in  the  absence  of  an  air  density 
channel  will  be  considered  first,  then,  the  effect  of  the  air  density  channel  on  the  projectile’s  net 
transverse  force  will  be  estimated. 

The  projectile  shape  is  a  two-dimensional  wedge  as  shown  in  Fig.  2.  The  characteristic 
transverse  length  is  Ja^  and  wedge  half-angle  is  0 .  The  projectile  is  assumed  to  be  traveling  at 
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supersonic  velocities  V  >  where  Vs  is  the  sound  speed.  An  oblique  shock  wave  is  produced 

having  an  angle  ft  with  respect  to  the  z  axis.  The  shock  angle  /i  is  a  function  of  the  Mach  number 
M  and  the  wedge  half-angle  0  .  The  shock  angle  is  always  greater  than  the  Mach  angle  // ,  i.e., 
J3  >  //  =  sin^fl/M),  where  the  Mach  number  is  given  by  M  =  V  /  Vs  (Anderson,  1991).  For  a 

stable  projectile,  a  body  downstream  from  the  wedge  is  required.  The  body  is  neglected  because 
it  would  have  complicated  the  analysis  of  the  projectile’s  net  forces  and  detract  from  the  focus  of 
this  paper,  which  is  on  the  characteristic  forces  from  an  air  density  channel  on  a  model  projectile. 


FIG.  2.  A  two-dimensional  wedge  projectile  shown  in  the  frame  moving  with  the  projectile.  The 
Mach  cone  (dotted  black  hnes),  fluid  (air)  streamlines  (red  arrows),  and  shock  wave  (solid  black 
lines)  are  also  shown. 

The  air  pressure  downstream  of  the  shock  wave  is  greater  than  the  ambient  air  pressure  and 
is  responsible  for  the  applied  normal  forces  on  the  projectile’s  surfaces.  By  summing  the 
components  of  these  forces,  the  drag  and  transverse  force  on  the  projectile  can  be  estimated. 

The  drag  force  is  defined  as  FD  =  (Cn  /  2)pA  V); ,  where  CD  is  the  drag  coefficient,  p  is 

the  mass  density  of  air  (p  =  1.2 kg/m3),  and  A±  is  the  characteristic  transverse  area  of  the 
projectile.  The  drag  coefficient  is  a  unitless  parameter  which  captures  any  additional  dependence 
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the  drag  force  may  have  on  Mach  number  M  ,  Reynolds  number  R  >  103 ,  and  projectile  geometry. 
Typically,  for  a  supersonic  projectile,  the  drag  coefficient  is  CD  ~  0.3  to  0.4  (McCoy  1990). 

If  the  ambient  air  density  is  uniform,  i.e.,  no  air  density  channel,  the  upper  (shown  in  Fig. 
2)  transverse  force  is  equal  in  magnitude  and  opposite  in  direction  of  the  lower  transverse  force, 
and  there  is  no  net  transverse  force  on  the  projectile.  There  would  be  a  drag  force  and  to  estimate 
it  we  consider  the  projectile  shown  in  Fig.  2.  The  oblique  shock  relations  (Anderson,  1991)  in  the 
limit  M2  sin2  j3  »l,  yield  a  pressure  on  the  upper  and  lower  surfaces  of  the  projectile 

pu  -  pt  -  2y(y  +  \ylpM2  sin“  J3 ,  where  p  is  the  upstream  pressure  and  y  =  1.4  is  the  heat 
capacity  ratio  of  air.  The  contributions  to  the  drag  force  from  the  upper  and  lower  surfaces  are 
Fd„  =Fd/  =-puAsin0z,  where  A  is  the  upper  or  lower  surface  area.  The  net  drag  force  is 

Fd  =  (CD  /  2)p  A  V2 ,  where  the  drag  coefficient  is  CD  -  4 (y  +  I)  1  sin2  /( ,  the  transverse  area  is 
A1=2Asin6l,  and  the  speed  of  sound  is  related  to  the  upstream  pressure  and  density  via 
Vs2  =  yp  /  p  ■  For  air,  CD  - 1 .7  sin2  j3 . 

If,  however,  the  ambient  air  mass  density  p  has  a  transverse  spatial  gradient,  such  as  an 
air  density  channel,  then  the  projectile  will  experience  drag  and  a  transverse  force.  Unforunately 
this  configuration  cannot  as  easily  be  analyzed  with  oblique  shock  theory  and  therefore  we  turn  to 
the  Newtonian  theory  is  used  to  gain  intutition  about  the  effect  of  an  air  density  channel.  The 
Newtonian  theory  neglects  shock  formation  and  calculates  the  nessecary  rate  of  change  in 
momentum  to  move  the  air  out  of  the  projectile’s  path.  For  a  wedge  in  uniform  air  the  drag  can 
be  found  by  using  the  oblique  shock  relations  and  taking  the  limits  of  M  — »  go  and  y—>l.  Asa 
result,  the  shock  angle  limits  to  the  wedge  angle,  /?—»(?  ,  and  the  drag  coefficient  becomes 
CD  =  2sin2  0  (Anderson  1991).  In  the  case  of  a  density  depression  of  A p  over  a  transverse  area 
A A  on  the  lower  side  of  the  wedge,  the  Newtonian  theory  yields  a  net  transverse  force  of 
h\  =  cot  0  CD  ( Ap  /  p)(AA/  A  ) p  A  V2  in  the  direction  of  the  air  density  channel,  where 

CD  =2  sin2  6k  In  order  to  get  an  estimate  for  the  characteristic  transverse  force,  we  can  use  a 
typical  drag  coefficients  of  CD=  0.3  (McCoy  1990)  to  provide  an  estimate  for  the  effective 
projectile  wedge  half-angle  6  ~  23° .  To  lowest  order  in  A p  and  A4 ,  the  ratio  of  the  transverse 
force  to  the  drag  is  F±  /  FD  =cot0(Ap/  p)(AA/  A±),  where  cot  6~  2.4  to  be  consistent  with 
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0~  23°.  While  the  fractional  change  in  air  density  and  ratio  of  air  density  channel’s  area  to 
projectile’s  area  are  both  small,  the  drag  force  can  be  sizable.  This  will  result  in  significant 
deflection  of  the  projectile. 

B.  Transverse  Deflection  of  Projectile 

The  net  transverse  force  on  the  projectile  will  cause  deflection  towards  the  air  density 
channel.  By  assuming  that  the  net  transverse  force  is  constant  as  the  projectile  passes  through  the 
air  density  channel  the  deflection  can  be  estimated.  A  single  air  density  channel  will  not  cause  a 
significant  deflection  by  itself.  Therefore,  several  air  density  channels  will  need  to  be  created  one 
after  another  to  build  up  the  transverse  velocity  of  the  projectile. 

For  a  single  air  density  channel,  the  transverse  velocity  of  the  projectile  v±(t)  is  given  by 

/ n  v±  =  F  ,  where  mp  is  the  projectile  mass  and  F  is  the  net  transverse  force  due  to  the  air  density 
channel.  As  will  be  shown,  the  initial  laser  pulses  can  be  changed  from  shot-to-shot  to  create  an 
extended  air  density  channel  that  is  formed  from  Nc  single  channels  joined  end-to-end.  The 
change  in  the  transverse  velocity  Av± ,  at  the  end  of  a  sequence  of  Nc  air  density  channels,  each 
of  length  L. ,  is  Av±  =  ( F  jmp  )( /V.L/V7^ ) .  Note  that  in  general,  the  channels  do  not  have  to 

connected  end-to-end  but  could  have  an  arbitrary  spacing. 

After  propagating  a  distance  L  from  the  start  of  the  air  density  channel,  the  transverse 
displacement  of  the  projectile  is  Ax±  =  Av±(L - NcLc  12)1  V p  ,  where  L>NcL  .  When  the 

estimate  of  the  net  transverse  force  is  used,  the  deflection  of  the  projectile  is  given  in  terms  of 
channel  and  projectile  parameters 

^  =  CjAgYMY  MxA'AV^A'Aj.  (1) 

2IpA4A  mr  A  2  ' 

As  an  example  of  the  amount  of  transverse  displacement  a  projectile  can  develop  we  consider  the 
parameters;  A/?//?  =  0.1,  A4/A±  =0.01,  CD  =0.3,  mp  =100g,  A±  =4cm2,  Nc  =20,  L  =  40m 

,  and  L  =  5km .  The  resulting  transverse  displacement  is  Ar  =2.9m.  For  comparison,  the 
deflection  after  just  the  first  air  density  channel  in  that  sequence  would  be 
Av1=(CD/2)(A/9/(,)(AA/Ai)(/>A14/m/,)(4/2)  =  0.6mm. 
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III.  Air  Density  Channel  Formation 

Formation  of  the  air  density  channels  occurs  in  a  four  step  process  (Cheng  2013,  Jhajj 
2014).  First,  an  ultrashort  laser  pulse  deposits  energy  in  the  air  via  multiphoton  ionization. 
Second,  the  free  electrons  redistribute  the  energy  gained  when  ionized  as  heat  by  colliding  with 
neutrals  and  attaching  to  oxygen  molecules.  The  local  heating  of  the  air  occurs  much  faster 
(nanoseconds)  than  the  energy  can  be  carried  away  by  acoustic  waves  (microseconds).  As  a  result, 
the  third  step  is  the  formation  of  a  pressure  shock  wave.  Finally,  after  the  region  has  reached 
uniform  pressure,  the  remaining  thermal  energy  has  an  associated  air  density  channel  which 
thermally  diffuses  on  millisecond  timescales. 

The  deposition  of  laser  energy  in  the  air  occurs  over  distances  larger  than  a  Rayleigh  length 
because  of  laser  filamentation.  Filamentation  is  the  result  of  a  dynamic  balancing  of  the  focusing 
effect  of  the  nonlinear  Kerr  index  and  the  defocusing  effect  of  the  laser  generated  plasma 
(Courairon  2007).  Due  to  this  balancing,  the  laser  intensity  and  plasma  density  in  an  air  filament 
are  clamped  to  approximately  1013 -1014  W/cm2  and  1— 5x  1016cm  3 ,  respectively  (Couairon 
2007).  Typical  electron  energies  in  a  filament  are  1  — 2eV  (Couairon  2007).  The  lifetime  of  the 
plasma  due  to  attachment  is  approximately  10ns  .  This  translates  into  plasma  filaments  which  at 
any  moment  are  approximately  3m  in  length.  These  filamentation  parameters  are  typical  and  have 
been  observed  both  experimentally  and  numerically.  They  will  be  used  when  making  estimates 
that  depend  on  the  filament’s  properties. 

The  details  of  the  laser’s  filamentation  (Sprangle  2002)  and  air’s  hydrodynamic  expansion 
(Cheng  2013,  Jhajj  2014)  have  been  discussed  in  previous  works.  Of  key  importance  to  projectile 
guiding,  are  estimates  of  the  magnitude  (AAO.),  volume  (  kR;L  ),  and  duration  (rdiff )  of  the  air 

channel’s  density  depression  and  the  resulting  requirements  placed  on  the  laser  pulses. 

To  estimate  the  magnitude  of  the  air  density  depression,  first  note  that  for  an  ideal  gas  at 
constant  pressure  the  change  in  air  number  density  AAC.  is  proportional  to  the  local  change  in  air 

temperature  ATT. ,  that  is,  AAC.  =  -AC.A77.  / 77. ,  where  the  total  number  density  and  temperature 

are  given  by  AC.  +  AAC.  and  77.  +  A 77. .  Estimating  the  increase  in  air  temperature  due  to  electron 

heating  will  provide  an  estimate  for  the  magnitude  of  the  air  density  depression.  For  a  gas  species 
in  thermodynamic  equilibrium  each  accessible  degree  of  freedom  (translational,  vibrational  and 


6 


NRL  Memorandum 


rotational)  has  a  kinetic  or  potential  energy  of  klsTa[  /  2 ,  where  kB  is  the  Boltzmann  constant. 
Therefore,  each  electron  and  air  molecule  have  a  thermal  energy  of  3 kBTe  /  2  and  5 kBTaii  /  2 , 
respectively.  Suppose  that  all  of  the  electron’s  thermal  (internal)  energy  is  transferred  to  the  air. 
The  change  in  air  temperature  due  to  electron  heating  is  (5  /  2)N,[u.kBATur  =  AUmr ,  where  the  energy 

deposited  in  the  air  is  A£/air  =  (3  /  2)NkBTe .  The  fractional  change  in  air  temperature  is 
ATair  /  rair  =  (3  /  S)(Ne  /  AC  )T  /  77.  *  0.085  for  7>2eV,  77,  =  0.026eV,  Ne=  5xl0]6cm'3,  and 
AC.  =  2.7xl019cm“3.  The  resulting  fractional  change  of  the  air  density  in  a  channel  is 
AAC  /  AC  -0.085. 

To  generated  air  density  channels  with  a  sufficient  density  depression  AAC.  for  projectile 
guiding,  enough  laser  energy  will  need  to  be  deposited  locally  to  heat  the  air.  This  sets  a  lower 
bound  on  the  necessary  laser  pulse  energy  £jasermin.  As  previously  mentioned,  the  temperature 

change  in  a  gas  is  related  to  the  thermal  (internal)  energy  by  A E  =  (5 /  2)NgjlkBATak  and  any 
changes  in  air  temperature  will  be  accompanied  by  a  change  in  air  density  AAC.  /  AC.  =  -A77r  / 77  r 
.  Therefore,  the  thermal  energy  change  in  a  single  air  density  channel  is  proportional  to  the  air 
density  depression  such  that  AE  =  (5/2)ACrkB77r(-AAC.  / ACr)(7z7?;L .) ,  where  R  and  L,  are  the 
channel  radius  and  length.  Typical  channel  dimensions,  as  will  be  discussed,  are  Lc  =40m  and 
Rc  =0.02 cm.  In  order  to  achieve  a  10%  depression  in  the  air  density  ( AAA.  / Nair  =-0.1),  the 
thermal  energy  must  increase  by  AE  =  140 mJ  .  This  sets  a  lower  bound  on  the  required  laser  pulse 
energy  Ejasermin  >  140mJ  and  is  achievable  with  current  laser  systems  (Stelmaszczyk  2004). 

The  air  density  channel’s  length  will  always  be  less  than  the  maximum  propagation  length 
of  the  filament  Lmax .  By  using  the  energy  depletion  length  of  the  filament,  an  upper  bound  on  the 
air  density  channels  length  can  be  estimated.  This  and  the  characteristic  channel  radius  of 
Rc  -  0.02cm  (Cheng  2013)  provide  an  estimate  for  the  volume  of  the  air  density  channel  kR] Lc  . 

The  maximum  propagation  distance  of  a  laser  pulse  generating  a  filament  is  limited  by  energy 
depletion  from  ionization  of  the  air  (Sprangle  2002).  If  all  of  the  laser  pulse  energy  goes  to  ionizing 
the  gas,  this  distance  is  approximately  Lmax  =  £jaser  /  {UvmN enR) ) .  For  a  laser  energy  of 

£jaser  =  140mJ  and  typical  atmospheric  filamentation  parameters,  such  as  an  ionization  energy  of 
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Ulon  =  12eV ,  electron  density  of  Ne  =  5  x  1016cm 3 ,  and  filament  radius  of  Rf  =  0.01cm  (Couairon 
2007),  the  maximum  propagation  distance  of  a  single  laser  pulse  undergoing  filamentation  is 
4a*  =46m. 

Finally,  there  are  several  relevant  timescales  to  consider.  The  first  is  the  timescale  for 
collisionally  heating  of  the  air  channel  by  the  electrons  (Tzortzakis  2001,  Cheng  2013).  The 
transfer  of  electron  kinetic  energy  to  the  neutral  gas  occurs  via  electron-neutral  collisions  and  the 
collision  time  is  rc  ~  lps .  As  a  result,  the  air  channel  heating  time  rh  is  approximately  the 
electron-to-neutral  mass  ratio  times  the  electron-neutral  collision  time  r. ,  that  is, 
rh  =  (M n  /  me)rc  ~  50ns ,  where  the  neutral  mass  is  Mn  ~  5x  10”26kg  .  This  is  comparable  to  the 

time  scale  of  electron  attachment  to  oxygen  and  it  is  expected  that  the  heating  mechanism  depends 
on  both  processes.  The  air  density  depression  develops  on  the  time  scale  that  the  pressure 
equilibrates,  that  is,  rdev  ~  Rc  /Vs  =0.1 /js  ,  where  the  speed  of  sound  is  Vs  =300 m/s  and 

Rc  =  200/mi .  The  air  channel  life  time  is  set  by  the  thermal  diffusion  rdiff  ~  R1  / 1)  =  \  ,5ms ,  where 

the  diffusion  coefficient  of  air  is  D  -  27mm2  /  s  (Cheng  2013)  and  R  =  200//m . 

As  an  aside,  the  electron’s  mean  free  path  is  small  compared  to  the  radius  of  the  filament 
implying  that  the  filament  and  air  channel  have  approximately  the  same  radius.  This  is  consistent 
with  experimental  observations  (Cheng  2013,  Courairon  2007). 

IV.  Flying  Focus 

The  concept  of  projectile  guiding  using  air  channels  relies  on  the  generation  of  many  air 
density  channels  along  the  projectile’s  trajectory.  As  discussed  previously,  the  characteristic 
lifetime  of  the  air  density  channel  is  a  millisecond.  For  a  projectile  velocity  of  Mach  3,  the  air 
channels  must  be  formed  no  farther  than  several  meters  in  front  of  the  projectile  or  else  the  air 
density  channel  will  have  dissipated  before  the  projectile  arrives.  Consequently,  the  plasma 
filaments  and  the  resulting  air  density  channels  must  move  from  shot-to-shot  at  the  projectile 
velocity. 

As  a  generalization  of  the  terminology  of  Gaussian  optics,  the  initial  location  at  which  the 
laser  pulse  collapses  into  a  plasma  filament  will  be  referred  to  as  the  focus.  The  movement  of  the 
focus  from  shot-to-shot  will  similarly  be  referred  to  as  the  “flying  focus.” 
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Linear  and  nonlinear  laser  pulse  propagation  provides  numerous  mechanisms  for 
controlling  the  laser  pulse  focus  or  filamentation  location,  such  as  focusing,  frequency  chirp,  or 
pulse  energy.  By  controlling  the  initial  laser  pulse  parameters  from  shot-to-shot,  the  focal  point 
of  the  laser  can  be  moved  and  thereby  create  a  flying  focus.  There  has  been  ample  theoretical  and 
experimental  work  on  ultrashort  laser  pulse  propagation  and  the  generation  of  filaments  in  the 
atmosphere  (Sprangle  2002,  Couairon  2007,  Berge  2007).  This  section  follows  the  methods  of 
previous  work  (Sprangle  2002)  in  solving  for  the  propagation  of  an  ultrashort  laser  pulse  in  the 
atmosphere. 

The  flying  focus  concept  considered  here  relies  on  propagation  of  intense  laser  beams  in  the 
atmosphere.  This  propagation  is  strongly  affected  by  various  interrelated  linear  and  nonlinear 
processes  (Sprangle  2002).  These  include  diffraction,  Kerr  self-focusing,  group  velocity 
dispersion,  spectral  broadening,  self-phase  modulation,  and  turbulence.  In  general,  a  laser  pulse 
propagating  in  air  can  be  focused  longitudinally  and  transversely  at  remote  distances  (~  km)  to 
reach  high  intensities  (~1013  W/cm2).  Due  to  group  velocity  dispersion,  temporal  pulse 
compression  can  be  achieved  by  introducing  a  frequency  chirp  on  the  pulse.  Nonlinear  transverse 
self-focusing  is  caused  by  the  optical  Kerr  effect.  Essential  to  filamentation  is  that  the  transverse 
self-focusing  leads  to  a  collapse  in  the  laser  spot  size. 

A  model  describing  the  longitudinal  and  transverse  compression  of  a  chirped  laser  pulses 
in  air  will  be  presented  (Sprangle  2002).  This  model  uses  an  extension  of  the  paraxial  wave 
equation  with  the  assumption  that  the  laser  pulse  maintains  a  Gaussian  shape,  both  transversely 
and  temporally,  as  it  propagates.  From  these  assumptions,  a  set  of  coupled  differential  equations 
for  the  laser  spot  size  and  temporal  duration  have  been  derived.  Additional  details  can  be  found 
in  reference  (Sprangle  2002). 

The  laser  electric  field  is  given  by  E(r,r/,r)  =  E(r,rj,z)  exp\iAkrj -icor]  ev /2  +  c.c. , 

where  E  is  the  complex  amplitude,  co  is  the  frequency,  k  is  the  wavenumber,  A k  —  k  —  cotc  is 
the  wavenumber  shift  from  air’s  linear  dispersion  relation,  r  is  the  radial  coordinate,  and 
t  -t-z/ c  and  /)  =  z  arc  the  light  frame  coordinates.  The  propagation  distance  z  and  time  t  are 
in  the  laboratory  frame.  The  complex  field  amplitude  is  assumed  to  be  a  transverse  and  temporal 
Gaussian  E(r,rj,r)  =  Ao  (77)  exp  [*#(77)]  expj^-(l  +  /a(^))r2  /  R2(r/)^  exp^-(l  +  z'/?(77))r2  /T2(t 7)J 

where  A0(tj)  is  the  field  amplitude,  0(ij)  is  the  phase,  R(r/)  is  the  spot  size,  a(i])  is  related  to 
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the  curvature  of  the  wavefront,  T(r/)  is  the  laser  pulse  duration,  and  [Mt])  is  the  chirp  parameter. 
The  quantities  Aq  ,  0 ,  T  ,  R  ,  a  ,  and  P  are  real  functions  of  the  propagation  distance.  As  shown 


in  reference  (Sprangle  2002),  the  following  coupled  equations  for  R  and  T  are  obtained, 


82R 


dr/2  k2R3 


1- 


K_}_ 

PnlTj 


(2a) 


&T  =  4P2E0  1  ,4  Pl 

dr/2 


k  PNL  RlT2 


(2b) 


where  E0  =  P(0)T (0)  is  proportional  to  the  laser  pulse  energy  and  is  independent  of  // , 
P(r/)  =  7tR2(j/)I(j/)  /  2  is  the  laser  power,  7(77)  =  cA;  (77)/ (8^r)  =  7(0)7?2(0)r(0)/[i?2(^)r(^)]  is 


the  peak  intensity,  PNL  =  A2  /  (2jinK)  is  the  self-focusing  or  critical  power,  and  /7,  is  the  group 
velocity  dispersion.  For  air  at  STP  and  wavelength  X  =  2n  I  k  « l//m,  the  group  velocity 
dispersion  is  /?2  =  2.2xl0“31s2 /cm,  the  Kerr  nonlinear  index  is  nK  =3x10  l9cnr  /  W ,  and 
1  +  nKI  is  the  refractive  index  of  air.  The  instantaneous  frequency  shift  along  the  pulse,  i.e.,  chirp, 
is  Sa>(r/,T)  =  2fi(r/)T /T2(j/) ,  where  fi(r/)  =  T(t/) / (2fi2)dT{r/) / dr/ .  A  negative  (positive) 
frequency  chirp,  /((//)  <  0  (/?(//)  >  0) ,  results  in  decreasing  (increasing)  frequencies  towards  the 
back  of  the  pulse.  Similarly,  the  wavefront  curvature  is  related  to  cc(j/)  =  —(kR(r/)  /  2)dR(r/)  /  dr/ 
with  a  focusing  (defocusing)  laser  field  corresponding  to  oc(r/)>  0  ( ait/)  <  0 ).  The  initial 
conditions  of  Eqs.  (2)  and  subsequent  evolution  of  R  and  T  are  determined  by  the  initial  laser 
spot  size  R( 0),  focusing  or  defocusing  via  a( 0),  pulse  duration  7T0) ,  and  frequency  chirp  via 


P(  0). 


The  first  term  on  the  right  hand  side  of  Eq.  (2a)  describes  vacuum  diffraction  while  the 
second  term  describes  nonlinear  self-focusing,  i.e.,  due  to  nK  .  Nonlinear  self-focusing  dominates 

diffraction  resulting  in  filamentation  when  P  >  PNL  ~3GW  (Sprangle  2002),  Couairon  2007). 

In  general,  the  laser  pulse’s  spatiotemporal  shape  will  not  remain  Gaussian,  but  this 
approximation  is  acceptable  during  the  propagation  leading  up  to  filamentation  (see  Fig.  10  of 
reference  [Sprangle  2002]).  After  filamentation,  the  laser  pulse  may  or  may  not  be  Gaussian  and 
ionization  would  need  to  be  included  in  the  model,  but  this  paper  is  only  interested  in  the  location 
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where  the  plasma  filament  begins.  Therefore,  solving  Eqs.  (2)  is  sufficient  for  determining 
feasibility  of  the  flying  focus. 

In  the  limit  that  the  pulse  length  does  not  change  appreciably  the  laser  spot  size  is  given  by 
R(r/)  =  R(0)[l-2a(0)T]/ZRO  +  (a2(0)  +  l-P/PA,i)(^/ZRO)2]  ,  where  ZR0=kR2( 0)/2  is  the 
Rayleigh  length  in  the  absence  of  initial-  or  self-focusing,  that  is  a(0)  =  0  and  PNL  — »  oo .  This  is 
similar  to  Gaussian  optics.  The  spot  size  collapses  to  zero,  R(r/c)  =  0,  at  a  distance 

rjc  /  ZRQ  -  (a(0)  -  yjp/ PNL  - 1  j  /  (a2  (0)  + 1  -  P  /  PNL  j  under  the  conditions  that  either 

(\+a2(0))PNI  >P>PNL  for  «(0)  >  0  or  P  >  (]  +  a2(0))PNI  for  all  a( 0) . 

A  similar  analysis  can  treat  the  laser  pulse  duration  without  self-focusing  ( PNL  — »oo) 

TOl)  =  T  (0)^(1  +  J3(0)tj  /Zt)2+(tj/  Zt  )2  ,  where  the  chirp  parameter  /?(())  is  like  the  initial 
temporal  focusing  or  defocusing,  and  the  group  velocity  length  ZT  =  7’2(0)/(2|/?2|j  is  like  the 
temporal  Rayleigh  length. 

For  example,  the  start  of  the  plasma  filament  is  shown  as  a  function  of  laser  pulse  energy 
in  Fig.  3.  In  the  ideal  case,  when  T(rj)  =  T( 0),  the  collapse  point  of  the  laser  spot  goes  off  to 
infinity  as  the  pulse  power  approaches  the  critical  power,  i.e.  0.6  mJ  for  a  200  fs  pulse.  Equations 
(2)  do  reproduce  the  result  of  z  =  ZR0  /  ^P  /  PNl  - 1  for  P  >  PNL  when  a(0)  =  0 .  but  when  the 

temporal  spreading  of  the  pulse  is  included,  higher  pulse  energies  are  requires  for  the  pulse  to 
collapse  at  the  same  propagation  distance.  With  temporal  spreading,  the  minimum  energy  for 
filamentation  is  approximately  three  times  larger.  The  collapse  point  of  the  laser  pulse  can  be 
moved  approximately  0.5  km  by  changing  the  pulse  energy  several  millijoules. 
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FIG  3.  The  filamentation  location  (or  collapses  point)  ijc  as  a  function  of  laser  pulse  energy  for 
T (0)  =  200  fs ,  R{ 0)  =  2  cm ,  a(0)  =  0 ,  J3(0)  =  0 ,  and  A  =  l//m .  The  red  dots  show  the 
numerical  solutions  of  Eqs.  (2)  and  the  red  solid  line  is  the  analytic  expression  for  r/c  when 
T  =  T( 0)  and  P  >  PNL .  The  black  dots  show  the  numerical  solutions  of  Eqs.  (2)  when  T ( z )  is 
allowed  to  evolve. 

The  flying  focus  concept  essentially  relies  on  changing  the  initial  conditions  of  the  laser 
pulse,  either  the  initial  energy  E0 ,  focusing  a(0) ,  chirp  parameter  /((()) ,  pulse  duration  T( 0) ,  or 
spot  size  R( 0) ,  from  one  laser  shot  to  the  next.  How  fast  the  flying  focus  can  be  moved  will 
depend  on  how  sensitive  the  filamentation  point  is  to  the  initial  conditions  and  on  how  quickly  the 
initial  conditions  can  be  changed  in  the  laser  system.  For  example,  if  the  laser  pulse  energy  is 
changed  shot-to-shot  and  we  neglect  GVD  (T(z)~  T0),  the  flying  focus  velocity  would  be  given 

by  vc  =  (drjc  /  dE0)(dE0  /  dt). 

For  projectile  guiding,  the  desired  flying  focus  velocity  is  ~1  km/s.  For  a  100Hz  laser 
system,  the  filaments  collapse  point  will  have  to  move  10  m  per  shot.  This  could  be  achieved  by 
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operating  a  laser  system  near  E0  ~  6  mJ  for  the  parameters  used  in  Fig.  3.  Near  this  energy  the 
rate  of  change  of  the  filamentation  point  is  drjc  /  dE0  «  45  m/  mJ .  The  laser  system  would  need  to 
be  able  to  change  the  pulse  energy  by  0.2  mJ  /shot. 

One  mechanism  for  experimentally  controlling  the  laser  pulse  energy  is  using  a  half-wave 
plate  and  thin  film  polarizer.  The  main  objective  of  a  system  like  this  is  to  take  the  laser  pulse 
after  the  last  laser  amplifier  and  dump  the  energy  that  is  not  needed.  For  example,  assume  that  the 
maximum  laser  pulse  energy  that  the  system  can  achieve  is  E0  max .  The  half-wave  plate  is  a 

birefringent  crystal.  If  the  crystal’s  length  is  properly  chosen,  then  the  relative  phase  delay 
between  the  components  of  the  field  along  the  crystal’s  optical  axes  will  cause  the  laser  field 
polarization  to  be  rotated  by  2 9  ,  where  0  is  the  angle  between  the  initial  laser  polarization  and 
the  crystal’s  fast  axis.  A  thin-film  polarizer  is  then  used  to  separate  out  one  polarization  and  create 
a  linearly  polarized  pulse  with  a  smaller  desired  energy  relative  to  the  maximum  that  the  system 
can  achieve,  that  is  E0  =  cos2  0EQ  max . 

For  example,  to  achieve  a  change  in  laser  pulse  energy  of  0.2  mJ/shot,  the  half-wave  plate 
can  be  rotated  at  the  rate  of  5.5  RPM  near  the  angle  6  ~  n  1 3.  The  laser  system  would  need  to 
have  peak  pulse  energies  of  E0  max  =  20  m.I .  to  create  5  mJ  pulses  in  this  configuration.  This  is 

approximately  what  would  be  needed  to  create  a  flying  focus  at  Mach  3  near  5  mJ  . 

V.  Conclusions 

In  conclusion,  we  have  detailed  the  mechanisms  for  guiding  a  supersonic  projectile  using 
optically  generated  air  density  channels.  An  estimate  of  the  transverse  force  and  deflection 
experienced  by  the  projectile  passing  through  air  density  channel  is  given.  We  have  outlined  the 
physical  processes  that  are  important  to  the  formation  of  a  channel  of  reduced  air  density. 
Additionally,  a  method  for  creating  extended  air  density  channels  by  using  a  flying  focus  was 
proposed. 
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